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APSIDAL MOTION OF THE ECLIPSING BINARY AS CAMELOPARDALIS: DISCREPANCY RESOLVED 

K. PAVLOVSKI^'^, J. SOUTHWORTH^ AND V. KOLBAS^ 
Draft version January 19, 2013 

ABSTRACT 

We present a spectroscopic study of the eclipsing binary system AS Camelopardalis, the first such 
study based on phase-resolved CCD echelle spectra. Via a spectral disentangling analysis we measure 
the minimum masses of the stars to be MAsin^i = 3.213±O.OO7M0 and M^sia^i = 2.323±0.006 M©, 
their effective temperatures to be Toff (A) = 12 840 ± 120 K and Teff(B) = 10 580 ± 240 K, and their 
projected rotational velocities to be WAsin«A = 14.5 ± O.lkms^^ and VBsinie ^ 4.6 ± O.lkms"^. 
These projected rotational velocities appear to be much lower than the synchronous values. We 
show that measurements of the apsidal motion of the system suffer from a degeneracy between orbital 
eccentricity and apsidal motion rate. We use our spectroscopically-measured e = 0.164±0.001 to break 
this degeneracy and measure Wobs = 0.133 ± 0.010° yr~^. Subtracting the relativistic contribution of 
= 0.0963 ± 0.0002° yr-i yields the contribution due to tidal torques: Wci = 0.037 ± 0.010° yr^^. 
This value is much smaller than the rate predicted by stellar theory, 0.40-0.87° yr~^. We interpret 
this as a misalignment between the orbital axis of the close binary and the rotational axes of its 
component stars, which also explains their apparently low rotational velocities. The observed and 
predicted apsidal motion rates could be brought into agreement if the stars were rotating three times 
faster than synchronous about axes perpendicular to the orbital axis. Measurement of the Rossiter- 
McLaughlin effect can be used to confirm this interpretation. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: early-type — stars: rotation 



1. INTRODUCTION 

Tidal torques dominate the dynamical evolution of 
close binary systems. In order of expected timescales, 
tidal effects should firstly align the stellar rotation axes 
with the orbital axis, synchronize the stellar rotational 
to the orbital frequency, and then circularize the orbit 
(|Mazehl I2008D . These effects happen on timescales of 
millions to billions of years, depending sensitively on the 
characteristics of individual systems. On a much shorter 
timescale apsidal motion, the precession of an eccentric 
orbit in its own plane, can produce an observable rate of 
change in the longitude of periastron, uj — duj/dt. 

AS Camelopardalis is a detached eclipsing binary con- 
taining two late-B stars with an orbital period of 3.43 d 
and a disputed apsidal motion period of order 2000 yr. 
The apsidal period, U, is much longer than expected 
for the measured properties of AS Cam, leading in the 
past to concerns about our understan ding of stellar 
physics and even of general relativity (iMalonev et al.l 
[1989, 1991). The apsidal motion of AS Cam was dis- 
covered by Khaliullin & Kozyreva (1983) from eclipse 
timings, and found to have a period oi U = 2250 ± 
200 yr. IMalonev et~aLl ()1989l) obtained a si milar U = 
2400t^ i^°o yr. Howev er, |KrzesinskLeLalJ (|1990D and 
iBozkurt fc Degirmencil (|2007. ') found very different apsi- 
dal periods, 920 ± 470 yr and 740 ± 6yr respectively, by 
adopting lower eccentricity values of e « 0.10. These 
shorter apsidal periods are closer to theoretical expec- 
tation s, but the low e confli cts with oth er observa- 
tions (|Hilditchl [19721 [Malonev et al.lll991|) . In addi- 
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tion, a third body has bee n found orbiting AS Cam with 
a period of P3 = 2.2 yr (iKozvreva fc Klialiullin[ 119991: 
IBozkurt fc Degirmencil I2007D . albeit with a low statis- 
tical significance. 

The binaries AS Cam and DI Herculis have for 
many years been the two best-known systems with 
apsidal periods much lo nger than they should be 
(iGuinan fc MalonevI Il985[ ). In a breakthrough work 
lAlbrecht et al.l ( 2000) observed the Rossiter-McLaughlin 
effect in DIHer and showed that its orbital and rota- 
tional axes are misaligned. This lengthens the expected 
apsid al period, which is now within 10% of the observed 
value ([Claret et al.ll201Ci[ ). Here we present the first high- 
resolution time-resolved spectroscopy of AS Cam, from 
which we find that the projected rotational velocities of 
the stars are much lower than expected and that their ro- 
tational axes are likely misaligned with the orbital axis. 
Below we refer to the primary (hotter and more massive) 
star as star A, and to the secondary as starB. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained 31 high-resolution spectra of AS Cam over 
five nights in 2007 October, using the Nordic Optical 
Telescope and high-resolution Fibre-fed Echelle Spectro- 
graph. Wavelength scales were established from thorium- 
argon exposures taken regularly each night. We opted 
for medium resolving power {R = 48 000) by using fibre 
3 in bundle B. This gave complete spectra coverage in 
the interval 3640-7360 A at a reciprocal dispersion rang- 
ing from 0.023 A px~^ in the blue to 0.045 A px~^ in the 
red, at a resolution of approximately 3.5 px. An expo- 
sure time of 600 s was used for all spectra, resulting in 
continuum signal-to- noise (S/N) ratios of roughly 90-140 
per pixel. 

Five additional spectra were secured in service mode 
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Fig. 1. — Example spectra of AS Cam in the region of tlie He I 4388 A line, offset from unity by arbitrary amounts for display purposes. 
The top four spectra (green lines) are observed composite spectra, with the orbital phases of observation annotated. The lower two spectra 
are the disentangled spectra for star A (blue line) and starB (red line). The disentangled spectra have been recombined and overplotted 
on the observed spectra (black lines) to show the quality of the fits. 

on 2010/10/10, in this case in the high-resolution {R = 
67000) mode with fibre 4 in bundle B. The spectral for- 
mat was identical to the medium-resolution mode but 
the resolution was 2.4 px. Exposure times of 900 s gave 
S/N ratios of about 250 per pixel. 

Basic reduction of these data (bias subtraction, flat- 
fielding, scattered light correction, extraction of orders, 
and wavelength calibration) was performed with IRAF. 
Removal of the instrumental blaze function was not triv- 
ial because the broad Balmcr lines of B-stars can extend 
over entire echelle orders. In such cases we interpolated 
between blaze functions from adjacent orders, using a 
semi-manual approach and JAVA routines written by VK. 
This process had to be undertaken carefully, as spectral 
disentangling is very sensitive to the normalization of the 
observed spectra. 



3. SPECTRAL DISENTANGLING ANALYSIS 

We subejcted our observations to a spectral disentan- 
gling (spd) analysis to determine the masses and atmo- 
spheric proper ties of the comp o nents of AS Cam. The 
SPD approach iSimon fc SturmI (|1994[ ) allows the best- 
fitting (in a least-squares sense) orbital elements and sep- 
arated spectra of a binary star system to be obtained di- 
rectly from a set of observed composite spectra covering 
a range of orbital phases (Fig.[T]). The resulting disentan- 
gled spectra have a very high S/N ratio as they contain 
the total signal in the input spectra (in this case about 
440 for star A and 180 for star B), so are well suited to fur- 
ther analysis such as chemical abundance determination 
jPavlovski & Hensbergc 2005|; iPavlovski fc SouthworthI 
[2009; Pavlovski ct al. 2009). spd does not require tem- 
plate spectra and is not biased by any blending of the 
spectral lines of the two stars in the observed spec- 
tra, so is excellent for measu ring orbital elements (e.g. 
iSouthworth fc Clauseiil l2007f ). A review of the theo- 
retical and practical aspects of SPD can be found in 



IPavlovski fc Hensbergd (|2009( ). 

We disentangled the spectra of AS Cam in F ourier 
space using the FDbinar\0 code (jlhiic et al.ll2004[ ). con- 
centrating on 10-15 relatively narrow spectral regions 
(50-80 A wide) which contain only metallic lines. Both 
stars exhibit numerous and very sharp spectral lines, re- 
sulting in well-defined orbital elements. Our adopted el- 
ements (Table[T]) are the mean and standard error of the 
results from the individual spectral regions. An impor- 
tant result is that eccentricity e — 0.164 ± 0.001, which 
conclusively rules out those apsidal motion studies which 
found short apsidal pe riods and e ~ 0. 10 (Sect.[T]). The 
velocity amplitudes of iHilditchI ()1972bf ). obtained using 
photographic methods, are in reasonable agreement with 
our values. 

We now turn to the measurement of effective tempera- 
ture (Toff) from Balmer line profiles. H/3 and H7 profiles 
for the two stars were obtained by disentangling the ob- 
served spectra with the orbital elements fixed at the val- 
ues in Table[TJ The disentangled spectra are not normal- 
ized to the correct continuum because such information 
is not available in t he input spectra (un less some were 
taken during eclipse: fTamaio et al.ll20lTI ). We therefore 
renormalized the Balmer profiles to the continuum using 
the light ratios of the two stars obtai ned by modellin g 
the UBV light curves of AS Cam from Iffilditcfl ()1972aD . 
The Teff values were then obtained by fitting the pro- 
files with synthetic H/3 and H7 spectra ca l cula ted using 
the UCLSYN program ( Smallevet al.ll200l]: ISmith 1992). 
We fixed the surface gravities of the stars to known val- 
ues, log5A = 4.154 ± 0.013 and logge = 4.278 ± 0.014 
(Southworth et al., in preparation). The blue and red 
sides of each Balmer profile were fitted separately, re- 
sulting in four Tog measurements for each star. Our fi- 
nal values are the mean and standard deviations of these: 

3 |http : //sail . zpf . f er . hr/f dblnary/ 1 
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TABLE 1 

The orbital elements and atmospheric properties of AS Camelopardalis from spectral disentangling. 



Parameter 



Star A 



StarB 



Orbital period (d) 
Time of periastron pasage (HJD) 
Orbital eccentricity e 
Longitude of periastron ui {°) 
Velocity amplitude (kms~^) 
Mass ratio q 

Minimum mass Af sin'^ i ( Mq) 
Effective temperature T^g (K) 
Projected rotational velocity vsini (kms~^) 



3.430973 (fixed) 
2454399. 7521f ±0.0026 
0.164 iO.OQf 
241.45 ±0.50 
106.22 ±0.14 146.92 ±0.13 

0.723 ±0.002 
3.213 ±0.007 2.323 ±0.006 

12 840 ±120 10 580 ±240 

14.5±0.1 ^4.6±0.1 
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Fig. 2. — Profiles of the Mgll 4481 A fine for star A (offset by 
+0.3) and starB (normalized to unity). Synthetic spectra for the 
measured vsini values are shown with red solid lines. Blue dot- 
ted lines show the synthetic spectra broadened to the synchronous 
rotational velocities. 

Toff(A) 12 840± 120K and rcff(B) = 10 580± 240K. 

The projected rotational velocities of the components 
were measured using a set of isolated metal lines, mostly 
of Fell, Tin and Cm, in the 4400-5000 A region. The 
instrumental broadening was obtained for each echelle 
order from the thorium-argon spectra. Representative 
values are 6.5 ± 0.1 and 4.8 ± 0.1 km s^^ for the 2007 
and 2010 data, respectively. Each line was then fit- 
ted with UCLSYN synthetic spectra, yielding ?;a sin ia — 
14.5±0.1kms~i and vb sinie ^ 4.6±0.1kms"^ These 
are much lower than the synchronous velocities, which 
are 36.7 zb 0. 5kms~^ a.nd 27 .1 ± 0.5 kms~^, respectively. 
In contrast, iHilditchI (|1972b) found that the vsini val- 
ues were approximately synchronous. We attribute this 
discrepancy to the long exposure times (roughly 1 hour) 
needed to obtain photographic spectra of AS Cam, which 
will have caused smearing of the spectral lines due to or- 
bital motion. Fig. [2] shows the Mgii 4481 A line profiles 
for the two stars compared to synthetic spectra calcu- 
lated for the measured and for synchronous v sin i values. 

4. APSIDAL MOTION 



Published measurements of the apsidal period and 
eccentricity of AS Cam are conspicuous by their dis- 
agreemt. Below we cast the argument in terms of the 
apsidal motion rate in degrees per year, cjobs is the ob- 
served rate, wgr is the relativistic contribution and Ud is 
the classical (tidal) contribution which can be estimated 
from stellar th eory. 

fKhaliuUin fc Kozyreval ()1983[ ) discovered apsidal mo- 
tion in the AS Cam syst em, finding e = 0.1695 ± 0.0014 
and Wobs = 0.16° yr-i. iMalonev et"all (| 19890 adopted 
e = 0.17 and got cjobs = 0.150 ± 0.053° yr"!, noting that 
this was far smaller than the expected amount due to 
c:;gr = 0.085 ±0 002° y r-^ and Wci = 0.358±0.058° yr"!. 
iKrzesinski et al.l (|1990[ ) suggested that a lower eccentric- 
ity could at least partially solve this problem , obtaining 
e = 0.10 ±0.01 and Wohs = 0.391^;;^^° yr"!. IWolf et alJ 
(|1996D adopted e = 0.14 and found Upbs = 0.18 3 ± 
0.026°yr-\ whereas iKozvreva k KhahullinI (|1999f ) as- 
sumed e = 0.17 to obtain c:;nhs = 0.1 49 ± 0.015° yr'^ 
Finally, Bozk urt fc Degirmenc I (|2007f ) arrived at e = 
0.1018 ± 0.0006 and Wobs = 0.486 ± 0.004° yr"!. This 
cjobs is consistent with wgr + ioci, although its errorbar 
is certainly too small. 

From the above details it is clear that e and Wobs are 
highly correlated for AS Cam, whose apsidal period is 
much longer than its observational history. To illustrate 
this we collected all available times of minimum light for 
AS Cam and analysed the m using the jktapsmot code 
(jSouthworth et aDl2004al ). which imple ments the exact 
ephemeris-curve method of lLacvl ()1992[ ). Most of these 
minimum timings have no associated errorbar, so un- 
certainties were estimated to be 0.01 d for photographic 
results and 0.001 d for photoelectric/CCD timings. All 
errorbars were then scaled up by a factor of five dur- 
ing the fitting process to force xl = 1-0 for the best 
fit. Solutions were then made both with (124 measure- 
ments) and without (92 measurements) the photographic 
timings, and for eccentricities of from 0.08 to 0.20. The 
upper panel in Fig.|3]demonstrates the strong correlation 
between e and Wobs, whereas the lower panel shows qual- 
ity of the fit. Whilst the lowest xt is found for e « 0.115, 
all solutions in the range e = 0.11-0.20 are statistically 
acceptable. 

We have a crucial advantage over previous stud- 
ies, namely the very precise e = 0.164 ± 0.001 found 
from our SPD analysis (Table [11). We can therefore 
reject low-eccentricity soluti ons (jKrzesinski et al.) 119901 : 
IBozkurt fc Degirmencil 120071) with extreme confidence, 
and also break the degeneracy between e and Wobs- Re- 
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Fig. 3. — Plot of the apsidal motion periods found for AS Cam 
as a function of eccentricity (upper panel). Filled red and open 
blue circles show results respectively with and without including 
photographic observations. The lower panel shows x2i which is 
forced to equal 1.0 for the best fit in each case. 

running the JKTABSDIM analysis with e fixed to 0.164 
and including all times of minimum, we find Wobs — 
0.133 ± 0.010° yr"\ corresponding to f/ = 2700±250yr. 

From the orbital elements of AS Cam (Table[T]) we mea- 
sure the relativistic contribution to the apsidal motio n 
to be WGR = 0.0963 ± 0.0002° yr"! (e.g. IGimene2jlT98l . 
Subtracting this from the observed value gives the rate 
due to classical effects: Wci = 0.037 ± 0.010° yr"!. 

The expected Wci can be obtained from the internal 
structure constant s log k-i , using the equations given by 
iClaret fc Gimene3 ()1993[ ). In turn, logfc2 must be es- 
timated from theoretical stellar structure models, and 
depends on the detailed characteristics of the stars. We 
have compared AS Cam to the Granada model tabula- 
tions (|Claretlll99,l [l997l : IClaret fc Giiniie^ [19951 [iQQSh 
using the minimum masses from Tablefl] orbital inclina- 
tion and radii from Southworth et al. (in preparation), 
and Tcff values from Sect. [3] From Fig.|4]we find that the 
fractional metal abundance of the binary is in the range 
Z = 0.01-0.02, which is qualitatively in good agreement 
with the spectral line strengths, and that its age is about 
130 Myr. 

The internal structure constants for Z = 0.02 and age 
105 Myr are logfc2(A) = -1.55 and logfc2(B) = -1.50, 
whereas for Z = 0.01 and age 160 Myr they are —2.33 and 
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Fig. 4. — Positions of the stars in AS Cam in the mass-radius 
and mass— Teff diagrams, compared to theoretical predictions from 
the Granada models. Results are shown for fractional metal abun- 
dances by mass o{ Z = 0.02 (green solid lines) and Z = 0.01 (blue 
dashed lines) and for the best-matching ages ±20 Myr. The filled 
circles represent the points at which tabulated model predictions 
are available, and the lines show a quadratic interpolation between 
these points. 



— 1.46. logfc2(A) is rather sensitive to the evolutionary 
status of star A. These logfc2 values result in theoreti- 
cal ujci values of 0.87 and 0.40° yr"\ for Z = 0.02 and 
Z = 0.01 respectively. These are both well in excess of 
the 0.037° yr"-'^ which is actually observed. We therefore 
confirm the discrepant apsidal motion of AS Cam to a 
high level of significance. 

The third body orbiting AS Cam has a very small 
light-time amplitude a n d an o rbital period of 2.2 yr 
(iKozvreva fc KhaliullinI 119991 : iBozkurt fc Degirmencil 
I2007f ). which is orders of magnitude shorter than U 
so should not have a significant effect on the analysis 
above. We have checked this by calculating a peri- 
odogram of the residuals of t he apsidal-mot ion fit, 
using the Period04 package ()Lenz fc Bregeii [2004). 
We find a peak of 2a significance at P3 = 824 d. This 
period is consistent with previous studies, as expected 
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because most minimum timings are in common. Our 
attempts to fit a spectroscopic orbit to the residuals 
led to a wide variety of solutions, depending on the 
starting parameter values and which parameters were 
fitted. New apsidal motion solutions with these orbits 
subtracted are not significantly different to our baseline 
solution. We conclude that there may be a third body, 
but that the data in hand are insufhcient to confirm its 
existence or give its orbital parameters, and that in any 
case it does not affect our measured Wobs- 

5. PROGNOSIS 

We find that AS Cam has a classical apsidal motion 
rate which is an order of magnitude lower than theoreti- 
cally predicted. How can this be explained? 

5.1. Problems with stellar theory or gravity theory 

The predicted Wci relies on our understanding of stellar 
physics, which is certainly imperfect. But it is difficult 
to see how the rate could decrease by an order of magni- 
tude in order to match our observations. Modern stellar 
theory does a good job of explaining apsidal motion in 
the g reat majority of close binaries (Claret fc Gimenezl 
12010^ so AS Cam would have to be a special case. Simi- 
lar comments apply to relativistic apsidal motion, as the 
theory of general relativity is otherwise highly successful 
in a wide range of scientific disciplines. 

5.2. Anomalously slow rotation 

The V sin i values of the stars are slower than the syn- 
chronous values by factors of 2.5 (star A) and 6 (starB). 
This apparently slow rotation results in a smaller pre- 
dicted Wcij but only by ~10%. Thus slow rotation cannot 
explain the Wci discrepancy. 

Tidal effects are expected to move the rotation rates 
towards the synchronous values and e towards zero, so it 
is reasonable to ask why we se e such low y sin i values. 
The synchronization timescales (jZahnlllQTSh are 106 Myr 
for star A and 1300 Myr for starB. As the age of the sys- 
tem is 100-160 Myr, we don't expect the stars to have 
reached synchronisation yet. They would thus have had 
to form with very slow rotation rates, which is excep- 
tional but at leas t more plausib le for a close binary than 
for a single star (|Tohlinell2002[ ). The orbital circulariza- 
tion timcscale is much longer again, in agreement with 
the observed eccentricity. 

5.3. Spin-orbit misalignment 

AS Cam and DI Her are two well-known binaries with 
anomalously low Wobs values which challenge our under- 
standing of stellar physics. The problem with DIHer 
has recently been explained as resulting from a large 
misalignme nt between the orbital and rotational axes 
(jAlbrecht et al. 2009. ). This is a highly plausible explana- 
tion for AS Cam, partic ularly due to the low y sin i values 
we find for both stars. iMalonev et al.l ()199lD considered 
this possibility for AS Cam itself, but rejected it because 
the stars were then thought to be rotating synchronously 
(jHilditch 1972b). 

In the light of the considerably subsynchronous rota- 
tion which we find for AS Cam , it is i nteresting to recon- 
sider this hypothesis. iShakural (|1985D found that the line 
of apsides can undergo retrograde motion if the stellar ro- 
tational axes are not aligned to the orbital axis, with the 



largest effect when the axes are perpendicular. In the 
case of AS Cam, Shakura found that the observed and 
computed did values could be brought into agreement if 
the stars were rotating at three times the synchronous 
rate around axes perpendicular to the orbital axis. It 
follows that the stars' rotation axes should be tilted by 
82° and 87° with respect to the orbital axis to produce 
the observed wsini values, which is indeed very close to 
perpendicularity. 

Although the axial misalignment hypothesis is very 
persuasive, it does incur the question of why tidal effects 
have not aligned the axes. The timescale for axial align- 
ment isrnuch shorter than for rotational synchroniza- 
tion (|Hutill981l : rMazehll2008D . and therefore much shorter 
than the age of AS Cam. Possible reasons for a mis- 
alignmen t between spin and orbital axes have been dis- 
cussed bv Bon nell et al.l (1992). Axial misalignment may 
also explain the slow rotational veloci ties found for the 
eclipsi ng binary systems V615Persei (jSouthworth et al] 
l2004bD and the central s tar of the planetary nebula 
SuWt2 (lExter et al.ll20Tol) . 

5.4. Perturbations from a third body 

The times of minimum light of AS Cam suggest (but 
do not require) the presence of a third body causing a 
light-time effect (Sect.|4|). Extensive numerical calcula- 
tions have shown that a third star in an orbit almost 
perpendicular to the orbital plan e of a close binary can 
cause anomalous apsidal niotion fehodvkin et al.l [20041 : 
iFabrvckv fc Tremainel[2007t iBorkovits et al.ll2007l) . This 
option seems unlikely, but cannot be dismissed as yet. 

6. SUMMARY 

We have presented the first modern spectroscopic 
study of AS Cam, an eclipsing binary which shows 
anomalously slow apsidal motion. Through a spectral 
disentangling approach we have obtained high-precision 
measurements of the minimum masses of the stars, their 
TcffS and their projected rotational velocities. If the 
stellar rotational axes are aligned with the orbital axis, 
the stars are rotating much more slowly than the syn- 
chronous velocities. We have re-investigated the apsidal 
motion of the system and demonstrated the strong corre- 
lation between apsidal motion rate and eccentricity. We 
have used our precise measurement of e = 0.1 64 ±0.001 to 
break this degeneracy, finding Wobs = 0.133±0.010° yr~^, 
corresponding to U — 2700 ± 250 yr. The relativistic 
component of this is wqr — 0.0963 ± 0.0002° yr~^, so the 
tidal component is therefore Wci = 0.037 ± 0.010° yr~^. 
From theoretical stellar evolutionary models we predict 
a very different uJc\ in the range 0.40-0.87° yr~^. 

We find no reason to suspect problems with our under- 
standing of stellar physics or general relativity, primar- 
ily because the Wobs values for most other close binaries 
agree well with theoretical predictions. Invoking slow ro- 
tation only changes the predicted Wd by ~10% so does 
not solve the discrepancy. However, the low usini values 
suggest that the rotational axes of the stars are highly 
inclined with respect to the orbital axes. S hakura (1985[) 
found that the discrepant apsidal motion for AS Cam 
could be explained if the rotational axes were perpendic- 
ular to the orbital axis and that the stars were rotating 
three times faster than synchronously. We therefore in- 
terpret our results as evidence of axial misalignment in 
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the AS Ca m system. The same phenomenon was found 
for DIHer ()Albrecht et al.|[2009l) . using a different obser- 
vational approach. 

Observation of the Rossiter-McLaughlin effect in 
AS Camelopardalis would allow further constraints 
to be placed on i t s dynamica l characteristics (e.g. 
lAlbrecht et all [20Q7L l2009l l20H . as would measuring 
additional times of minimum light. Further theoretical 
study of the tidal effects in misaligned binary systems 
would also be very illuminating. 
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